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Abstract 

Parameter fitting of data to a proposed equation almost always consider these 
parameters as independent variables. Here, the method proposed optimizes an ar- 
bitrary number of variables by the minimization of a function of a single variable. 
Such a technique avoids problems associated with multiple minima and maxima 
because of the large number of parameters, and could increase the accuracy of the 
determination by cutting down on machine errors. An algorithm for this optimiza- 
tion scheme is provided and applied to the determination of the rate constant and 
final concentration parameters for a first order and second order chemical reaction. 



1 Introduction 

Deterministic laws of nature are sometimes written - for the simplest examples- in the 
form 

Yia^^Yia^{F,k,t) (1) 

linking the variable Yiaw to t. The components of P, Pi{i = l,2,...Np) and k 
are parameters. Verification of a law of form ([1]) relies on an experimental dataset 
{{Yexp{ti), ti),i ~ 1,2, ...N)}. Confirmation or verification of the law is based on 
(a) deriving suitable values for the parameters (P, k) and (b) showing a good enough 
degree of fit between the experimental set YexpiU) and Yiaw{ti)- Many methods 
|[l]|2][3]lll etc.] have been devised to determine the optimal P, k parameters, but most if 
not all these methods consider the aforementioned parameters as autonomous and inde- 
pendent (e.g. [2 1) subjected to free and independent variation during the optimization 
process. On the other hand, if one considers the interplay between the experimental 
data and Yiaw one can derive certain parameters like the final concentration terms (e.g. 
Aoo and Yoo in what follows in sec© ) if k, the rate constant is known. To preserve 
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the viewpoint of the inter-relationship between these parameters and the experimental 
data, we devise a scheme that relates P to fc for all Pi via the set {Yexp{ti),ti}, and 
optimize the fit over fc-space only. i.e. there is induced a Pi{k) dependency on k via the 
the experimental set {Yexp{ti), ti}. it is unclear at present whether this optimization 
procedure is equivalent to previous ones, but its structure is not in contradiction with 
situations where there are inter-relations between the variables, and the results for the 
first and second order kinetics presented here are in very close agreement with those 
derived from the published literature, the advantages of the present method is that the 
optimization is over 1 D k space, leading to a unique determination of P with respect 
to k, whereas if all P are considered equally free, the optimization could lead to many 
different local solutions for each of the {Pi}. In what follows here, we assume that the 
rate laws and rate constants are not slowly varying functions of the reactant or product 
concentrations, which has recently from simulation been shown generally not to be the 
case Ii5j. 



2 Outline of Method 

As above, N is the number of dataset pairs {YexpiU), ti}, Np the number of compo- 
nents of the P parameter, and Ng the number of singularities where the use od a partic- 
ular dataset (Yexp, i)leads to a singularity in the determination of Pi{k) as defined be- 
low and which must be excluded from being used in the determination of Pi{k). Then 
{Np + l) < {N~N,) for the unique determination of {P,fc}. Define ^-^'Cat^ = 
as the total number of combinations of the data-sets {Yexp{ti), ti} taken Np at a time 
that does not lead to singularities in P^. Write Yiaw in the form 

Yia^^it,k)^ f{P,t,k). (2) 

Then map / — > YthiP, t, k) as follows 

Yth{t,k)^ f{P,t,k) (3) 

where the term P and its components is defined below and where fc is a varying pa- 
rameter. For any of the (ii, 12, . . . , iat ) combinations where ij = {Yexp{ti .),ti .) is 
a particular dataset pair, it is in principle possible to solve for the components of P in 
terms of k through the following simultaneous equations: 

Yexp{ti;^) — /(P,iij,fc) 



Yexp{ti2) — f(P,ti2,k) (4) 

For each Pi, there will be Nc different solutions, Pi{k, 1), Pi{k, 2), . . . P,;(fc, iVJ . We 
can define (there are several possible mean definitions) an arithmetic mean for the 
components of P as 
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Each Pi{k,j) is a function of k whose derivative is known either analytically or by 
numerical differentiation. To derive an optimized set, then for the least squares method, 
define 

N' 

Q{k) = J2{Yexp{U)-Yth{k,U))\ (6) 
1=1' 

Then for an optimized k, we have Q'{k) — 0. Defining 

N' 

Pk{k) = - Yth{k,U)).Y;f^{k,U) (7) 

1=1' 

the optimized solution of k corresponds to Pk{k) — 0. The most stable numerical 
solution is gotten by the bisection method where a solution is assured if the initial 
values of k yield opposite signs for Pk{k). Since all Pi{k) functions are known, their 
values may all be computed for one optimized k value of Q in (|6]l. For a perfect fit 
of Yexp with Ylaw, Q{k') — Q'{k') = ^ Pj ^ Pj (Vj) and so in this sense we 
define the above algorithm as giving optimized values for all Pj parameters via the 
k determination. This method is illustrated for the determination of two parameters 
in chemical reaction rate studies, of l^*- and 2^^^^ order respectively using data from 
published literature , where this method yields values very close to those quoted in the 
literature. 

3 Applications in Chemical Kinetics 

The first order reaction studied here is 

(i) the methanolysis of ionized phenyl salicylate with data derived from the literature 
|lllTable7.1,p.381] 

and the second order reaction analyzed is 

(ii) the reaction between plutonium(VI) and iron(II) according to the data in ||7] Table 
II p. 1427] and lE Table 2-4, p.25]. 



3.1 First order results 

Reaction (i) above corresponds to 

PS^+CHgOH ^ MS"+PhOH (8) 

where the rate law is pseudo first-order expressed as 

rate = fcJPS]" = fcc[CH30H][PS-]. 

with the concentration of methanol held constant (80% v/v) and where the physical 
and thermodynamical conditions of the reaction appears in 16] Table 7.1,p.381]. The 
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change in time t for any material property X{t), which in this case is the Absorbance 
A{t) (i.e. A{t) = X{t)is given by 



A(t) ^\oo- (Aoo - Ao) exp i-kj) (9) 

for a first order reaction where Aq refers to the measurable property value at time t = 
and Aoo is the value at t = oo which is usually treated as a parameter to yield the 
best least squares fit even if its optimized value is less for monotonically increasing 
functions (for positive ^at all t) than an experimentally determined X{t) at time t. In 
Table 7.1 of |6| for instance, A(t = 2160s) = 0.897 > Aopt,oo = 0.882 and this value 
of A^ is used to derive the best estimate of the rate constant as 16.5±0.1 x 10~^sec~^. 
For this reaction, the Pi of (|2]i refers to Aoo so that P = Aoo with Np = 1 and k = ka- 
To determine the parameter Aoo as a function of ka according to (|6]l based on the entire 
experimental {{Xexp, U)} data set we invert (|9]l and write 

(fc) = ± y (Ae.pfe)-A,exp-fct.) 

where the summation is for all the values of the experimental dataset that does not 
lead to singularities, such as when ti — 0, so that here A^^ = 1. We define the non- 
optimized, continuously deformable theoretical curve Af ^ where \th = , fc) in © 

as 

\th(t, k) = Aoo(fc) - (Aoo(fc) - Ao) exp (-M) (11) 

With such a projection of the Aoo parameter P onto fc, we seek the least square min- 
imum of Qi(fc), where Qi(fc) = Qof® for this first-order rate constant k in the 
form 

N 

Qlik) = Y.'^XexpiU) - XthiU, k) f (12) 

i=l 

where the summation is over all the experimental {Xexp{ti), ti) values. The resulting 
Pfc function Q for the first order reaction based on the published dataset is given in 
Fig. ©.The solution of the rate constant k corresponds to the zero value of the function, 
which exists for both orders. The P parameters (Aoo and Foo ) are derived by back sub- 
stitution into eqs. (fTOl l and ( fTSl l respectively. The Newton-Raphson (NR) numerical 
procedure |9, p.362]was used to find the roots to P^.For each dataset, there exists a 
value for Aoo and so the error expressed as a standard deviation may be computed. The 
tolerance in accuracy for the NR procedure was 1. x 10^^" . We define the function 
deviation fd as the standard deviation of the experimental results with the best fit curve 
fd = Vw{Y.iLii^exp{ti) - Xth{ti f} Our results are as follows: 
ka = 1.62 ± .09 X lO-^s-i; Aoo = 0.88665 ± .006; and fd = 3.697 x lO^^. 
The experimental estimates are : 

ka = 1.65 ± .01 X lO^^s-i. = 0.882 ± 0.0; and fd = 8.563 x 10"^^ 
The experimental method involves adjusting the Aoo = Aoo to minimize the fd func- 
tion and hence no estimate of the error in A^o could be made. It is clear that our method 
has a lower fd value and is thus a better fit, and the parameter values can be considered 
to coincide with the experimental estimates within experimental error. Fig.Qshows 
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the close fit between the curve due to our optimization procedure and experiment. The 
slight variation between the two curves may well be due to experimental uncertainties. 
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Figure 1 : Plot of the experimental and curve with optimized parameters showing the 
very close fit between the two. The slight difference between the two can probably be 
attributed to experimental errors. 



3.2 Second order results 

To further test our method, we also analyze the second order reaction 

Pu(VI) + 2Fe(II) Pu(IV) + 2Fe(III) (13) 

whose rate is given by rate = /co[Pu02^][Fe'^^] where fco is relative to the constancy 
of other ions in solution such as The equations are very different in form to the 
first-order expressions and serves to confirm the viability of the current method. 

For Espenson, the above stoichiometry is kinetically equivalent to the reaction 
scheme |^ eqn. (2-36)] 

PuO^++Fe2+ A PuO++Fe3+ 

which also follows from the work of Newton et al. [7, eqns. (8,9),p.l429] whose 
data El TABLE II,p.l427] we use and analyze to verify the principles presented here. 
Espenson had also used the same data as we have to derive the rate constant and other 
parameters [8 , pp. 25-26] which is used to check the accuracy of our methodology. The 
overall absorbance in this case Y{t) is given by |8 , eqn(2-35)] 

_ + {^0 (1 - a) - Foe} CXp(-fcAof) 

^'^ l-aexp(-fcAoO ^ ^ 
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Figure 2: functions ^ for reactions (i) and (ii) of order one and two in reaction 
rate. 

where a = is the ratio of initial concentrations where [B]o > [A]o and [B] = 

[Pu(VI)], [A] = [Fe(II)] and [B]o = 4.47 x lO^^M and [A]o = 3.82 x IQ-^M . A 
rearrangement of ( fT4b leads to the equivalent expression ^ eqn(2-34)] 

\ [A]o (Ft - yoo) J [A]o ' 

According to Espenson, one cannot use this equivalent form [8, p. 25] "because an 
experimental value of Yrx, was not reported." However, according to Espenson, if Yoo is 
determined autonomously, then k the rate constant may be determined. Thus, central to 
all conventional methods is the autonomous and independent status of both k and 1^00- 
We overcome this interpretation by defining Yoo as a function of the total experimental 
spectrum of ti values and k by inverting (fT4l i to define Yoo{k) where 

Y^{t)= ^ yexpfe){exp(fcAot,)-a)} + yo(a-l) ^^^^ 



N' (exp(fcAot,) - 1) 



where the summation is over all experimental values that does not lead to singularities 
such as at ti = 0. In this case, the P parameter is given by Yoo(fc) — Pi{k), fct = k is 
the varying k parameter of (|2]). We likewise define a continuously deforming function 
Yth of k as 

_ Yoojk) + {Yq (1 - g) - Yooik)}eM-kAot) 
1 — a exp(— fcAotj 
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In order to extract the parameters k and Yqo we minimize the square function Q2{k) 
for this second order rate constant with respect to k given as 

N 

Q2{k) = J2^Ye^P^U) - YthiU, k)f (18) 

i=l 

where the summation are over the experiment ti coordinates. Then the solution to the 
minimization problem is when the corresponding Pk function (I?) is zero. The NR 
method was used to solve Pk — with the error tolerance of 1.0 x 10^^". With the 
same notation as in the first order case, the second order results are: 
h = 938.0 ± 18M s-i; Foo = 0.0245 ± 0.003; and fd = 9.606 x lO^''. 

The experimental estimates are JH p.25]: 

h = 949.0 ± 22 X lO-^s-i; = 0.025 ± 0.003. 

Again the two results are in close agreement. The graph of the experimental curve and 
the one that derives from our optimization method in given in Fig.©. 
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Figure 3: Graph of the experimental and calculated curve based on the current induced 
parameter-dependent optimization method. 



4 Conclusions 

The results presented here show that for linked variables, it is possible to derive all 
the parameters associated with a curve by considering only one independent variable 
which serves as a function of all the other variables in the optimization process that 
uses experimental dataset as input variables in the estimation. Apart from possible 
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reduced errors in the computations, there might also be a more accurate way of deriving 
parameters that are more determined by the value of one parameter (such as k here) than 
others; the current methods that gives equal weight to all the variables might in some 
cases lead to results that would be considered "unphysical". 
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